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ABSTRACT 

During 2003 and 2004 the Anomalous X-Ray Pulsar XTE J1810-197 went through a series of 
four bursts. The spectrum in the tail of one of these bursts shows a strong, significant emission 
feature 13 keV, thereby encoding a wealth of information about the environment surrounding 
this object. In this paper we analyse this emission feature considering both cyclotron and atomic 
emission processes and weigh our findings against three leading AXP models: the Magnetar 
model, Fall-back disk model and the Quark nova model. We find that atomic emission from 
Rubidium within a Keplerian ring (~15 km from a compact object of ~ 2M Q ) is the most 
consistent scenario with the observations, supporting the Quark nova model. Cyclotron emission 
from an atmosphere a few hundred meters thick also fits the feature well, but is ruled out on 
account of its positional coincidence in three separate AXP sources. 

Subject headings: stars: magnetic fields - stars: neutron - stars: individual (XTE J1810-197) - X-rays: 
bursts 



1. Introduction 

Anomalous X-Ray Pulsars (AXPs) are solitary 
objects characterized by persistent emission, X- 
ray pulsations and a rapid spin down rate. Rarely, 
some AXPs also emit Soft Gamma-ray Repeater 
(SGR) -like X-ray b ursts. The first detected AXP 
was IE 2259-586 bylFahlman & Gregoryl (fl981|. 



Since then, several more have been discovered in- 
cluding IE 1048.1-5937, 4U 0142+61, IE 1547.0- 
5408, 1RXS J170849.0-400910, IE 1841-045, XTE 
J1810-197, CXOU J010043.1-721134, and CXO 
J164710.2-455216. These objects exhibit spin pe- 
riods of 5-12 seconds, fairly constant x-ray lumi- 
nosities of ~ 1 x 10 33 erg s _1 to ~ 1 x 10 36 



erg- 



soft x-ray spectra and rapid spin do wn rates 



of P = 1 x 10" 



-13 



1 X 10 



" 10 s s- 1 



(Woods 



et al. 2005). AXPs have some properties similar 



to SGRs, and are thought to be a related phe- 
nomenon. 

The leading model to explain the AXP phe- 



nomenon is the Magnetar model (Thompson & 
Duncan| [l995 ) . Magnetars are neutron stars (NS) 
with extreme magnetic fields between ~ 1 x 10 14 
and ~ 1 x 10 15 G. The bursting mechanism in 



the Magnetar model may be located on the NS 
surface, or in the magnetosphere. In one model 
the large magnetic field exerts stress on the crust, 
causing it to fracture and subsequently emit short 
X-ray bursts. These "starquakes" are believed to 
create additional hot-spots on the surface of the 
NS (Pons & Rea 2012). Another model suggests 



that reconnection events of a highly twisted mag- 
netic field create the bursts in the magnetosphere 
(Lyutikov et al. 2002). The spin-down rate and 



periods observed in AXPs are explained through 
a combination of magnetic-breaking and magnetic 
field decay on a time-scale shorter than the spin- 
down time-scale. 



Another AXP model was put forth by Chat 



terjee et al. (2000) in which a fossil accretion 
disk supplies the X-ray luminosity. The disk is 
formed by fall-back material from the supernova 
(SN) explosion that created the NS. The Fall-Back 
Disk (FBD) model re quires modest magnetic fi elds 
(B <~ 1 x 10 13 G, |Ertan fc Erkut] (|2008|; 



Er- 



tan et al. (2007)) with inner and outer radii of 
~10,000 km and >100,000 km respectively. It 
should be noted that the FBD model does not 



explain the origin of the bursts and some hybrid 



magnetar models have been suggested ( Ek§i & Al- 



par 2003 Ertan & Alpar 20031. Evidence of a 



fall-back disk in the source 4U 0142+61 has been 



reported by Wang et al. ( 2006 ) . Ertan & Caliskan 
[ ( 2006 ) provide further evidence of fall-back disks 
around other AXPs giving support to this the- 
ory. This spin-down rates and periods observed 
in AXPs are explained through torques provided 
by the magnetic-field/disk interaction. (Ertan et 
n|2009l ) 



al 



In recent years a new theory to explain the 
properties of SGRs and AXPs has emerged using 
quark stars as the engine ( [Ouyed et al~1|2007a|b| ) . 
A quark star is created through the detonation 
of the neutron star in a process called a quark 



nova (QN) (Ouyed et al. 


2002 


Keranen et al. 


2005 Niebergal et al. 


2010 


I. The QN ejects the 



neutron-rich outermost layers of the neutron star 



(Ouyed & Leahy 2009), leading to r-process pro- 



duction of heavy elements ( Jaikumar et al.[2007 |. 
The QN can result in heavy-element-rich mate- 
rial from the ejected NS crust in a Keplerian ring 
orbiting the remnant quark star, extending from 
~20 km to ~100 km depending on its age( Ouyed 



et al. 2007b). This ring is distinguished from a 
fall-back disk, as it is much closer to the surface 
of the central object, and is made up of primar- 
ily degenerate material. The quark star is born 
with a dipole magnetic field aligned with its axis 
of rotation ( Ouyed et al.|2004 2006 1 and a surfac e 
magnetic field of B ~ 1 x 10 Jijl G ( |Iwazaki||2005 1 . 
In this model, X-ray bursts from AXPs occur when 
chunks of material break off the ring and accrete 
onto the quark star surface, while quiescent emis- 
sion is powered by magnetic flux expulsion during 
spin-down. In this model, the period and spin- 
down rate of AXPs is through dipole spin-down 
and magnetic field decay (Niebergal et al. 2006). 



AXP Spectra 

During AXP quiescent phases, their spectra be- 
low ~10 keV are often featureless with a contin- 
uum characterized by the sum of a power law 
and black body, the sum of two black bodies or 



even by resonant cyclotron scattering (Reaetal 
2008[ ) . Broac[] emission lines have been observed 



1 Gavriil et"al] ( 2008| repo rt "very broad" line s in all three 
sources in their analysis; |Woods et al.| ( |2005| > report a 3<r 



in the spectra of three AXPs during their burst- 
ing pha se; 4U 0142+61, XT E J1810-197 and IE 
1048.1 ( |Gavriil et al.||2008j ). In each case, there 
is a significant feature at ~13 keV which has the 
potential to provide crucial information about the 
bursting environment. In this paper we will con- 
centrate on the spectrum from XTE J1810-197. 

XTE J1810-197 w as first observed in 2003 by 
Ibrahim et al. (2004 1. It is at a distance of ^5 



kpc ( Gotthelf et al.||2004 ), has a period of 5.5 sec- 
onds and a spin down rate of P ~ 1 x 10~ n . 
The inferred magnetic field due to magnetic brak- 
ing is - 3 x 10 14 G. In 2005, X TE J1810-197 wen t 
through a series of four bursts (IWoods et al.|2005|) . 



During one of the burst tails, a significant emis- 
sion line at ~13 keV was detected. Subtraction 
of the best-fit continuum model leaves a strong, 
significant emission feature at ~13 keV (Figure 9 
from |Woods et al.| ( [2005])). |Gavriil et al.|(|2008| re- 
analysed the data of Woods et al. ( 2005[ ) intheir 
comparison of spectral features from 4U 0142+61, 
XTE J1810-197 and IE 1048.1. Their re-analysed 
data is shown in Fig[l] with a linear energy axis. 
This work will focus on exploring what physical 
and geometrical conditions are required to repro- 
duce the width and location of the 13 keV feature 
of XTE J1810-197 and the implications for AXP 
models. 

In the Magnetar model, the line is com- 
monly interpreted as a cyclotron line, implying 
field strengths ~ 2.2 x 10 15 G for protons and 
- 1.2 x 10 15 G for electrons ( |Gavriil et aL]|2008 |. 
In this work we argue that the large width of the 
13 keV feature is naturally explained (see section 
3.2 below) by cyclotron radiation in a spatially 
varying magnetic field. The width of the line then 
provides information on the extent of the emit- 
ting region in a dipole or multipolar field around 
the magnetar. The cyclotron interpretation does 
have some challenges: if due to protons, the in- 
ferred magnetic field (~ 2.2 x 10 15 G) is much 
larger than that derived by magnetic breaking 
(~ 3 x 10 14 G). This can be reconciled with non- 
dipolar fields however, because the spin down rate 
depends on the dipolar component, whereas the 



upper limit of 4.1 keV width for XTE J1810-197. The dif- 
ferent interpretations in line width is related to detector 
limitations, such as limited sensitivity and spectral resolu- 
tion, and to different data analysis methods. We proceed 
here using the |Gavriil et ah| | |2008[ ) (Fig. 2) spectral data. 
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line gives the field at the surface. If the line is an 
electron cyclotron line, the implied field is too low 
for the surface of the magnetar, and thus requires 
emission from far above the surface (~100 km due 
to the nature of a dipole field). Another com- 
plication is that there are other significant lines 
detected (at ~4 and -8 keV in 4U 0142+61). 
If these are also proton cyclotron lines, then the 
emission must be coming from multiple, distinct 
regions on the magnetar surface, possibly provid- 
ing evidence for a multipolar field. Perhaps the 
most damaging argument against the cyclotron 
interpretation of the 13 keV line is that it occurs 
in the three sources at the same energy, implying 
that these sources all have an emission region with 
the same magnetic field. 

Another possible interpretation is that the 13 
keV feature is an atomic line. This resolves the 
problem of why we see the same line in three dif- 
ferent sources. One of the challenges, as usual, is 
identifying which species is responsible for the line. 
The high energy of the feature implies a heavy 
(heavier than Fe), highly ionized atom, which may 
be difficult to explain in some models. Another 
challenge with the atomic line interpretation is 
the width of the feature. Doppler broadening re- 
quires velocities a good fraction of the speed of 
light which are not expected in or around a mag- 
netar with a period of 2-10 s. 

The purpose of this paper is to investigate the 
origin of the 13 keV feature by considering both 
cyclotron and atomic lines in different models. 
The paper is organized as follows: In section 2 
we discuss the modelling procedure and the soft- 
ware we use. Section 3 investigates the viability 
of the cyclotron interpretation by investigating the 
emission in several different scenarios. Section 4 
then examines the atomic line interpretation and 
its feasibility. Section 5 discusses the results and 
the implications for the different AXP models. Fi- 
nally Section 6 provides our conclusions and pre- 
dictions for future observations. 

2. Line Modelling with SHAPE 

In order to investigate both cyclotron and 
atomic emission from an astronomical source, 
we need software capable of performing radiative 
transfer. Additionally, because we are exploring 
such a vast parameter space, we need the software 



to be fast and easy to use. Software that meets 



these two requirements is SHAPE (Stcffcn et al 



2011 



SHAPE is a morpho-kinematic modelling tool 
developed to model astronomical objects and their 
environment in three dimensions (3D). The envi- 
ronment within SHAPE allows one to easily con- 
struct 3D structures with arbitrary density, tem- 
perature and velocity profiles. SHAPE then calcu- 
lates the resulting emission and displays output in 
the form of two dimensional (2D) images, PV di- 
agrams and spectra for comparison with observed 
data. 

In order to keep the processing time at a min- 
imum, SHAPE uses a ray-tracing radiative trans- 
fer algorithm. The difference between this ap- 
proach and other radiative transfer code is that 
SHAPE requires certain attributes (e.g. tempera- 
ture, density, velocity, emitting species etc..) be 
defined as input, instead of them being calculated 
by the code. SHAPE therefore answers the ques- 
tion "Given a temperature, density and velocity 
distribution, what is the resulting spectrum from 
a particular species?" 

It should be noted that SHAPE does not take 
radiative effects specific to high magnetic fields 
into account. For example, birefringence, photon 
splitting, scattering suppression and distortion of 
atoms in large magnetic fields are not taken into 
account. 

2.1. Algorithm and Spectral Modelling 

SHAPE partitions 3D space into an axis aligned 
(with the observer) grid of 256 3 cells. Each cell 
within the volume of a structure (e.g. sphere, 
torus, disk) is assigned attributes (e.g. temper- 
ature, density, velocity, pressure, magnetic field, 
gravitational field) from which emission and ab- 
sorption coefficients are calculated depending on 
the species involved (e.g. atoms, molecules, elec- 
trons, protons). A ray is then passed through 
the grid for each pixel of the detector (computer 
screen) travelling through multiple cells where ra- 
diation is either added to, or removed from the 
ray. Each ray consists of m frequency bands evenly 
spaced between a minimum and maximum fre- 
quency range, set to match the observed data. 



http: / /www.astrosen.unam.mx/shape/ 



3 



The energy within each band of each ray can then 
be plotted, generating a final spectrum that can 
be compared with observations. For more infor- 
mation on the algorithms involved in SHAPE, see 
http ://www.astrosen.unam.mx/shape/! 

For the cyclotron simulations, SHAPE assigns 
each cell in the grid a magnetic field based on the 
selected field type (dipole or constant) and its lo- 
cation in space. The field within each cell is there- 
fore taken as constant and we do not consider the 
effects of a charged particle travelling from one cell 
to another. 

For each of the models presented in this work, 
spectra will be produced with SHAPE in the afore- 
mentioned manner. The energy range will span ^3 
- 30 keV in order to match the observed spectra 
of XTE J1810-197 and divided into 100 frequency 
bands. When comparing our best-fit models with 
the data, we apply a Gaussian convolution with 
FWHM of 1.5 keV to simulate the energy resolu- 
tion of the observations. A grid of 256 3 will be 
used in all models. 

3. Cyclotron Emission 

Cyclotron radiation is the result of a moving 
charged particle being deflected in a static mag- 
netic field. Due to the Lorentz force, a moving 
charged particle will undergo circular motion per- 
pendicular to its velocity and the magnetic field. 
For particles that are moving only slightly rela- 
tivistically, radiation is primarily emitted at the 
cyclotron frequency: 



eB/jm 



(1) 



Where e is the charge of the particle, B is the 
magnitude of the magnetic field, 7 is the relativis- 
tic correction and m is the mass of the charged 
particle. Lines can also be emitted at harmonics 
of the fundamental frequency, their line centres 
given by = I x cj c , where I is an integer rang- 
ing from I — > 00. The line center of the emission 
is therefore proportional to the magnitude of the 
magnetic field. The emission coefficients for the 
cyclotron lines are given by |Boyd fc Sanderson~| 
p)03l: 



ne 2 l 2 uj 2 l 21 
167r 2 e c[Z!] 2 



(cos^ + ljsin 2 * 1 - 1 ^) 



knT P 



2mc 2 



2nk B T e 



1/2 



exp 



—rnc 2 {u — uji) 2 
2k B T e uj 2 cos 2 (9) 



(w, cos(fl))- 1 
(2) 



Where n is the density of the charged parti- 
cles, I is the harmonic number, 9 is the angle be- 
tween the magnetic field and the observer, k B is 
the Boltzmann constant, T e is the electron temper- 
ature, and m is the mass of the charged particle. 
In the regime where k B T e /2mc 2 -C 1, the funda- 
mental frequency will dominate the spectrum. For 
this reason we will concentrate on the fundamen- 
tal line, 1 = 1. The emission coefficient then takes 
the simpler form: 



ne 2 u 2 c k B T e /cos 2 (6>) + 1 



32m7r 2 e c 2 \ lo c cos{9) 
—mc 2 (u! — lo c ) 2 
2k B T e uj 2 cos 2 (6») 



2nk B T e 



exp 



Equation [3] represents a Gaussian line profile 
with amplitude: 



ne 2 uj 2 k B T e ( cos 2 (6>) + 1 



32m7r 2 eoc 2 \ u; c cos(#) 
and width: 



2irk B T e 



1/2 



Aw FWHM 



2.35fc B T e w 2 cos 2 (9) 



-1I/2 



(4) 



(5) 



For a given temperature and density, the inten- 
sity of the line is proportional to the angle of the 
observer and the magnetic field, and the mass of 
the charged particle: 



a oc 



B 

^5/2 



>s 2 (6>) + l 
cos(0) 



(6) 



This implies that as the angle increases, so will 
the intensity and a proton is expected to produce 
a less intense line than an electron. The width of 
the line is also proportional to the angle between 
the observer and the magnetic field and the mass 
of the charged particle: 



Aw 



FWHT 



Bcos(9) 

m 5/2 



(7) 



1/2 



(3) 
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Smaller angles will create a broader line, and 
electrons will produce broader lines than protons. 
For a given line center, uj c , we can expect the elec- 
tron line to be ~ 43 x broader than that produced 
by protons. 

We wish to investigate the cyclotron line shape 
for both electrons and protons in a variety of en- 
vironments in an effort to resolve the origin of the 
13 keV emission feature seen in XTE J1810-197. 

3.1. Constant Magnetic Field 

We begin by investigating the trivial case of 
thermalized electrons and protons in a constant 
magnetic field. The geometry of the model is that 
of a spherical shell with an inner and outer radius 
of 10 km and 12 km respectively. The density 
within the shell is constant, as is the temperature 
which is set to 1 x 10 8 K. 

We model both proton and electron cyclotron 
radiation with a constant magnetic field of 2.15 x 
10 15 G and 1.15 x 10 15 G respectively. We use 
these field strengths so that the line centres match 
up with the 13 keV feature in the spectrum of 
XTE J1810-197. Since the line shape depends on 
the angle between the magnetic field and the ob- 
server, we will examine two scenarios; one with 
the magnetic field at 30° to the observer and one 
at 60°. 

The left panel of Figj2] shows both the electron 
and proton results. The most dramatic result be- 
tween the proton and electron lines is their relative 
breadth. The breadth in the proton case, at both 
30° and 60°, does not exceed one frequency band 
in our simulation. On the other hand, the electron 
cyclotron lines are much broader. This is of course 
due to the mass dependency of the line width (see 
equation [5]). The electron line under the 30° incli- 
nation is more intense but narrower than the 60° 
line. This again can be predicted by equations [4] 
and [5] The width of a cyclotron line is therefore 
heavily dependent on both the mass of the parti- 
cle and on the angle between the observer and the 
magnetic field. 

From this simulation we can see that the 
breadth of the 13 keV feature can easily be ob- 
tained through electron cyclotron emission in a 
constant magnetic field. For proton cyclotron 
emission, the larger mass of the proton severely 
limits the line breadth and another mechanism is 



required for its broadening. 

3.2. Constant Magnetic Field + Variation 

A perfectly uniform, constant magnetic field is 
not realistic and some variation should be ex- 
pected. The next simulation investigates the 
shape of the proton cyclotron line in a non- 
uniform, variable magnetic field. For this test, the 
magnitude of the magnetic field at each location 
within the shell will be randomly chosen between 
2.15 x 10*5 G ±dB , and oriented at 60° to the 
observer. We will run two simulations, one with 
dB = 0.5 x 10 15 G and one with dB = 1.0 x 10 15 
G. 

The right panel of Figj2]shows the results of this 
second simulation. As expected, as the variation 
in the magnetic field (dB) increases, the breadth 
of the line increases. A variation in the magnetic 
field of dB = 1.0 x 10 15 G is sufficient at 60° for 
proton cyclotron radiation to reproduce the width 
of the 13 keV feature. 

3.3. Dipole Magnetic Field 

The previous simulations deal with the ideal- 
ized case of charged particles in a (relatively) con- 
stant magnetic field. In realistic situations, such as 
the space around a magnetar, the magnetic field 
configurations are more complex. It is tempting 
to conclude that because the cyclotron line center 
is proportional to the magnetic field, a continu- 
ous, variable magnetic field like that of a mag- 
netic dipole will produce a flat continuum of ra- 
diation, and not a distinct line. The situation is 
complicated, however, by the dependence of the 
line strength on the angle between the field and 
the observer. The next set of simulations is meant 
to investigate this situation by studying the cy- 
clotron line properties produced from a spherical 
shell within a magnetic dipole field. 

The parameters of this simulation are the same 
as before; a spherical shell with an inner and 
outer radius of 10 km and 12 km respectively and 
a constant temperature of 1 x 10 8 K. The an- 
gle between the magnetic moment of the dipole 
and the observer is fixed at 90°. For the pro- 
ton cyclotron case, the magnetic dipole moment 
is /i m = 2.0 x 10 30 A m 2 such that the magnetic 
field inside the spherical shell is ~ 3 x 10 15 G and 
within our spectral range. For the electron cy- 
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clotron case, we use /x m = 1.0 x 10 27 Am 2 . To 
show the dependency of the line shape on the size 
of the emission region, we have plotted the emis- 
sion from a shell 1 km and 2 km in thickness, with 
the outer radius in both cases held at 12 km. 

The left panel of Fig|3] shows the results for the 
electron case and the right panel the results for 
the proton case. The dipole field introduces sev- 
eral levels of complexity over the constant mag- 
netic field case which is represented by the spectra. 
Each point within the spherical shell has a mag- 
netic field with a different direction and magni- 
tude. Since the cyclotron lines are highly sensitive 
to the angle between the field and the observer, 
the location and intensity of the lines from each 
location in space can be highly variable. We can 
see this variability most clearly in the proton cy- 
clotron spectrum in the form of "spiky'j^] features. 
Interestingly, however, the sum of each individual 
line contributes to a broad, asymmetrical feature 
and not a flat continuum. The electron case, as 
expected from the greater width of each individ- 
ual line, shows a broader feature with less vari- 
ability. In both the electron and proton cases, the 
thinner shell model resolves two distinct features, 
whereas the thicker shell blends the features into 
one broad one. These broad features, especially 
when binned to instrument resolutions, may ex- 
plain one or more features in the AXP spectra. 

3.4. Dipole Magnetic Field Polar Regions 

It is entirely likely that the emission of the cy- 
clotron lines does not come from a complete spher- 
ical shell, but a small region perhaps around the 
poles of a dipole or multipole field. This next sim- 
ulation examines the consequence of the cyclotron 
lines emitted in a region of space near the poles 
of the dipole field. We use the same configuration 
as before, except the spherical shell is limited to 
a region between = 0° and <j) = 20° where <fi 
is the angle measured from the pole (</> = 0°) to 
the equator (0 = 90°). As with the previous sim- 
ulation, we simulate both a thicker shell (2 km) 
and a thinner shell (1 km). The results are shown 
in the left and right panels of Fig|4] for the elec- 
tron and proton cases respectively. We can clearly 
see that in both cases, the thickness of the shell 



3 The spectrum shows these spiky features due to the finite 
size of the individual grid cells in our simulation 



contributes to the thickness of the emission fea- 
ture. The width of the observed 13 keV feature can 
therefore be matched in both the electron and pro- 
ton cases with emission from a thin shell around 
the polar region of a dipolar magnetic field. 

3.5. Dipole Magnetic Field Polar Regions 
- Rotated 

The 13 keV feature is fit quite well by either 
proton or cyclotron emission from a small region 
around the poles of a magnetic dipole field. We are 
aware, from our previous simulations, of the effect 
of field orientation on the intensity and width of 
the line. This next simulation therefore builds on 
the previous by examining the consequences of ro- 
tating the dipole field with respect to the observer. 
We run four simulations, with the magnetic mo- 
ment at 90°, 60°, 30° and 0° degrees to the ob- 
server. Fig[5]show the results for both the electron 
and proton cases. As the angle between the dipole 
moment and the observer decreases, the line ap- 
pears to shift to the red. Interestingly it does not 
get much broader, but roughly keeps its original 
shape. This may seem counter intuitive since the 
position of the line, according to equation [T] de- 
pends only on the magnitude of the magnetic field. 
This is true in a constant field, but in this situation 
we have a range of fields so we have emission over 
many frequencies. It is the orientation that deter- 
mines which lines are enhanced and suppressed. 
At steeper angles, the polar regions are enhanced 
(which also have the strongest magnetic field and 
highest energy lines). As the angles decrease, the 
polar regions are suppressed and the surrounding 
regions, which have lower magnetic fields, are en- 
hanced resulting in an emission feature at lower 
energies. 

3.6. Conclusion 

The above simulations were meant to show how 
the cyclotron emission line looks under various 
conditions. We have shown that the width and 
intensity of the line depends heavily on both the 
mass of the charged particle and the angle be- 
tween the magnetic field and the observer. Further 
adding to the breadth of the line is the variation 
in the magnetic field. Looking at the realistic case 
of a polar emission region in a dipolar field, we see 
that the breadth of the line depends more on the 
thickness of the emission region than the mass of 



G 



the particle. Somewhat surprisingly, the direction 
of the magnetic dipole moment with respect to the 
observer along with its magnitude determines the 
center of the emission feature. 

We can now put all of this information together 
and try to fit the 13 keV feature seen in XTE 
J1810-197. Figj6] shows the best fit models for 
the electron and proton cases in the left and right 
panels respectively. In both cases we have a dipole 
magnetic field and the emitting region is a spheri- 
cal shell restricted to the angles <fi = 0° — > 20° . In 
both cases, we require a thin shell with an inner 
and outer radius of 11.9 km and 12 km respec- 
tively. For the electron case, the magnetic dipole 
moment is \i m — 1.0 x 10 27 A m 2 inclined at 66° 
to the observer. For the proton case, the magnetic 
dipole moment is /i m = 2.0 x 10 30 A m 2 , inclined 
at 56° to the observer. 

Of course these fits are not unique as the same 
line profile can be achieved in several ways. For 
example, the line center can be shifted by chang- 
ing the angle of the dipole field with respect to the 
observer, or simply adjusting the magnitude of the 
magnetic dipole moment. What these fits do tell 
us is that both electron and proton cyclotron emis- 
sion can explain the 13 keV feature seen in XTE 
J1810-197 under reasonable physical and geomet- 
rical conditions. 

4. Atomic Line Emission 

In this section we will investigate the idea that 
the 13 keV feature in the spectrum of XTE J1810- 
197 may be caused by atomic line emission. There 
are several questions concerning the atomic line in- 
terpretation. First, is there a viable atomic tran- 
sition near 13 keV and second, how do we explain 
the large breadth of the 13 keV feature. By fit- 
ting the 13 keV feature with a Gaussian profile, 
we find that the center of the line is at 13.07 keV, 
the standard deviation is a sa 1.45 keV, and the 
full-width-half-max is AE 1/2 w 3.41 keV. 

We use the National Institute of Standards and 
Technology (NIST) atomic line database to iden- 
tify possible counterpart lines to the 13 keV fea- 
ture of XTE J1810-197. Tabled] lists all the avail- 
able transitions between ~9 keV and ~18 keV 
with determined transition probabilities. 

The most likely species for the counterpart line 
is Rubidium (Rb). We come to this conclusion 



for several reasons. First, and most important, is 
that only Rb lines have sufficient energy to be a 
viable candidate. There are four Rb transitions 
near 13 keV: 14.3 keV (RbXXXVII), 14.2 keV 
(RbXXXVII), 13.9 keV (RbXXVI) and 13.8 keV 
(RbXXVI). The next closest candidate is Ni with 
an energy of 10.356 keV, which is much too low. 
Note that all the Rb lines are actually above 13.07 
keV. This is not a problem since when these lines 
are gravitationally red-shifted, they will fall closer 
to the observed feature. Second, the Rb lines at 
~13.8 keV all have the highest transition proba- 
bilities of any line in the entire spectrum (lines 
spanning the whole spectrum were not listed in 
the table due to the large number of transitions) . 
This and the fact that several of these lines are ex- 
tremely close, makes a strong feature at ~13 keV 
very likely. 

Fig[7] shows the line identifications of the tran- 
sitions listed in Table [T] plotted with the observed 
spectrum from XTE J1810-197. It is clear from 
this figure that Rb is the most likely candidate for 
the 13 keV feature. Interestingly several weaker 
Rb transitions lie close to the energy of a small 
"bump" around 17 keV. This bump, however, may 
or may not be real. 

Two ionization states of Rb may contribute to 
the observed feature; RbXXXVI and RbXXXVII. 
The ionization energy of RbXXXV (giving RbXXXVI) 
is 4.36 keV an d that of RbXXXV I (giving RbXXXVII) 

Using E 



is 18.31 keV (Sansonetti 



2006). Using E = kT 
temperatures on 
2.0 x 10 8 



we find that this translates to 
the order of T w 5.0 x 10 7 K and T 
K for RbXXXVI and RbXXXVII respectively. 
For this study we will assume a temperature of 
T = 1.5 x 10 s K. We are unsure of the electron 
density, so we assume an abundance of RbXXXVI 
10% that of RbXXXVII. 

4.1. Line Modelling 

We will model the emission from Rb assuming 
it is optically thin and that we are in LTE. The 
emission coefficients for a bound-bound transition 
from level i — > j are calculated via: 



3v 



4-7T 



(8) 



Where is the number density of atoms in 
the ith state and A^j is the transition probability. 
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In LTE we can use the Boltzmann equation to 
determine the number density for state i given a 
temperature T and total density n. Employing 
equation [8j for the Rb transitions in [T] we can use 
SHAPE to produce a model spectrum. The left 
panel of FigjH] shows the resulting spectrum for a 
temperature of T = 1.5 x 10 8 K and a resolution 
of 1000 frequency bands in order to capture the 
finest structure. 

4.2. Gravitational Red-shift 

With the most likely species selected, we can 
now go ahead and investigate the physical and 
geometric conditions required for the Rb line to 
reproduce the 13 keV feature seen in XTE J1810- 
197. We will use the Rb XXXVII transition at 
14.299 keV as our reference line (most intense line 
in the left panel of FigjHJ . 

The gravitational red-shift near a massive ob- 
ject is given by: 



E Q = E eX l- 



2GM 



(9) 



Where E a is the observed energy, E e is the 
emitted energy, G is the gravitational constant, 
M is the mass of the compact object and r is the 
distance from the object. Solving for r, we can 
determine the location of the emission above the 
compact object: 



2GM 



El 



ei - m 



(10) 



For a mass of 1.4M Q , and a shift from 14.299 — ► 
13.07 keV, we get an emitting radius of 25 km from 
the center of the compact object. If we instead 
take the closer line, Rb XXXVI at 13.794 keV, we 
get an emitting radius of 60 km. If the 13 keV 
feature in XTE J1810-197 is indeed an Rb atomic 
line, it must therefore be emitted from at least 25 
km from the center of a compact object of 1.4M© 
(assuming relatively small velocities so relativistic 
beaming is not a factor). 

We model this scenario by placing a spherical 
shell with an inner and outer radius of 24.5 km 
and 25.5 km respectively in a gravitational field 
produced by a compact object at the origin. The 
right panel of FigjH] shows the resulting spectrum, 
the strongest line matching perfectly with the cen- 
ter of the 13 keV feature. 



4.3. Line Broadening 

Next we will look at what conditions can lead to 
an atomic line as broad as that seen in XTE J 1810- 
197. The two most likely broadening mechanisms 
for the 13 keV feature are Doppler broadening due 
to the bulk motion of the atoms, and thermal 
broadening due to their random motion. Using 
the relativistic Doppler formulation, the line-of- 
sight velocity corresponding to a shift in energy 
E e - E is: 



{El -El) 
'{El+E*) 



At the 2cr level, E a — 15.97 keV, we have veloci- 
ties of V2<j ~ 0.20c and at the 3cr level the velocities 
reach as high as v^ a 0.29c. If the width of the 
line is due to non-relativistic, Maxwellian thermal 
broadening, this translates into a temperature of 
QKwok |[2#7] : 



c 2 a 2 
k R E^ 



(12) 



The lowest possible temperature (using hydro- 
gen for the mass) is therefore T « 1.1 x 10 11 
K. For Rubidium we require a temperature of 
T « 9.5 x 10 12 K. These temperatures are ex- 
treme and not expected in any AXP model. We 
will therefore not consider this option further. 

Since these speeds are a significant fraction of 
the speed of light, relativistic beaming is expected 
to be important and may alter the line shape. The 
intensity at a specific frequency is scaled as: 



In = D A E 



(13) 



Where I and I e are the observed and emitted 
intensities respectively,!? = 1/7(1— j3 cos(0)) is the 
Doppler factor, and 9 is the angle between the ob- 
server and the particle velocity. Due to relativis- 
tic beaming, we expect to see the blue emission 
increased, and the red suppressed assuming mo- 
tion toward/away from the observer. This should 
also cause the peak of the line to shift towards the 
blue. 

There is another possible broadening mecha- 
nism that can be significant near a massive, com- 
pact object. Near such an object, the gravita- 
tional field may significantly shift the line center 
towards the red as given by equation [9j If the 
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emission region is extended, however, the mate- 
rial closer to the mass will experience a greater 
red-shift than the material further away. The net 
effect will be "Gravitational Broadening" of the 
spectral line with the breadth conveying the ex- 
tent of the emission region. 

With the above considerations, we will now 
model these different broadening mechanisms in 
an effort to determine under what conditions we 
can get a broad enough line using atomic transi- 
tions. 

4--3.1. Doppler Broadening 

We will investigate Doppler broadening of the 
Rb line under two different geometries: a spheri- 
cal shell and extended disk. The velocity profile 
will be one of rotation, where the velocity direc- 
tion is perpendicular to the y-axis and the position 
vector (observer line-of-sight is along the z axis). 
The magnitude of the velocity will be constant at 
0.2c. The density is constant, as is the tempera- 
ture which is set to T = 1.5 x 10 8 K. We use 100 
frequency bands. 

Spherical Shell 

The first model we investigate is a spherical 
shell with an inner and outer radius of 24.5 km 
and 25.5 km respectively. The left panel of Figj9] 
shows the results of this simulation. We incline 
the sphere (rotational plane) at 90° and 30° to the 
observer to demonstrate the effects of inclination 
on the spectrum, with 0° being pole-on. As pre- 
dicted, the width of the lines roughly match that 
of the 13 keV feature and the Doppler broaden- 
ing suppresses the red while enhancing the blue. 
In the 90° rotation case, this leads to a shift of 
the line center towards the blue. The 30° rota- 
tion of the sphere means the line of sight velocity 
is decreased, and we therefore see less of a shift. 
Because of the relativistic beaming, the center of 
the line in the 90° rotation case is shifted towards 
higher energies. In order to match the 13 keV 
feature, we require a larger gravitational red-shift 
which means the emitting region must be closer 
than the 25 km calculated above. 



The next model we investigate is a rotating 
disk. The disk has an inner radius of 24.5 km, 
an outer radius of 25.5 km and a height of 5 km. 
We model the system at two inclinations: 90° and 
30° , with the results displayed in the right panel of 
Figj9] The resulting line is broadest at an inclina- 
tion of 90° and narrowest at 30° . This is expected 
since the line of sight velocity is maximum when 
the inclination of the disk is 90°. It is also ap- 
parent that the feature is being split at the lower 
inclinations. This is a common spectral profile of 
a rotating disk. The red-shifted line is lower than 
the blue-shifted because of relativistic beaming. 

4-3.2. Gravitational Broadening 

If emission is from an extended atmosphere in a 
large gravitational field, it is possible for a spectral 
line to experience gravitational broadening. Using 



equation 10 we can get an idea of the extent of the 
atmosphere. 

Immediately we run into problems with this 
analysis. The Rb line is emitted at 14.299 keV 
but the feature is observed to extend blue-ward 
of this. Gravitational red-shift only shifts the line 
to the red, making this scenario impossible (again 
assuming relativistic beaming is negligible). How- 
ever, there are weaker Rb lines at greater ener- 
gies (~ 16 — >-~ 19 keV) which may be contribut- 
ing to the extended breadth of the line. For this 
reason we will consider the maximum possible en- 
ergy of the line under gravitational broadening to 
be 14.299 keV. Under this assumption, the line 
therefore extends from E re( ^ ~ 11.84 — > E-^^q = 
14.299 keV. For a compact object of 1.4M this 
translates, through equation|10| to an inner radius 



of r-: 



13 km. The outer radius, in theory, must 



extend to infinity. If we take a more reasonable 
maximum energy of the line to be -Eblue ~ 14.2 
keV we find the outer radius must be at least 
7'i n ~ 300 km. 

We notice, by equation |9j that most of the en- 
ergy of the resulting feature will be emitted closer 



4 A "disk" has different meanings for each of the AXP mod- 
els. In the QN model, the disk is made up heavy-element- 
rich degenerate material and located ~20 km to ~100 km 
away from the quark star. The disk of the FBD model is 
formed from the debris of the supernova. It is much larger 
and located further away from the central star. 
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to the non-shifted energy (more emitting atoms 
are at larger radii) . For the red side of the feature 
to be observable, a density gradient is therefore 
required. We model this situation in SHAPE by 
using a disk with an inner and outer radius of 10 
km and 300 km respectively in a gravitational field 
of a 2Mq compact object (we use 2M instead of 
1.4M to make the effect more pronounced) We 
model three different density distributions, the re- 
sults of which are shown in FigjlOj Panels a, b 
and c show the line profile using n = constant, 
n oc 1/r 2 and n oc 1/r 3 respectively. We can see 
that as the density profile becomes steeper, the 
line becomes broader. This is because there is 
more of the inner material where the largest red- 



shift is occurring. Panel d of Fig 10 shows the re- 



sults of the n oc 1 jr A profile but convolved with a 
Gaussian with FWHM of 1.5 keV to better match 
the observations. We can see that the breadth 
and location of the line can be reproduced in this 
scenario. 

4.4. Conclusion 

The 13 keV feature of XTE J1810-197 may 
be explained by line emission from Rb XXXVI 
or RbXXXVII. If the feature is indeed from an 
atomic line, the atoms must be either moving 
at mildly relativistic speeds, emitted from an ex- 
tended atmosphere in a large gravitational field, 
or perhaps a combination of both. Although the 
feature is probably a blend of many Rb lines, the 
14.299 Rb XXXVII line is the most likely candi- 
date for line center. In order to explain the offset 
in energy between this Rb line and the observed 
feature, the emission must be coming from a re- 
gion at least 25 km from the origin assuming a 
compact object of 1.4M Q . We have seen, how- 
ever, that when large velocities (0.2c in this case) 
are involved, relativistic beaming can shift the line 
center blue-ward. To shift the line center back to 
the location of the observed feature, the gravita- 
tional red-shift must be enhanced either by moving 
the emission region closer in, or by increasing the 
mass of the compact object. 

We can take this knowledge and create a best 
fit to the observed feature of XTE J1810-197. For 
a rotating shell, we get a best fit with a spherical 
shell having an inner and outer radius of 20 km 
and 25 km respectively, a compact object mass of 
M = 1.4M© , a constant velocity of v = 0.2c and 



an inclination (of the rotation plane) of 32°. The 
results are shown in Fig|ll| Again this fit is not 
unique as, for example, a rotation in inclination 
is equivalent to decreasing the velocity. The fit 
reproduces the width of the 13 keV feature quite 
well and also reproduces the small "bump" at ~17 
keV. As mentioned earlier, this bump may or may 
not be real. In either case, the Rb atomic line 
model does predict it at the given temperature. 

For the rotating disk case, we found the opti- 
mal parameters to be a disk with an inner and 
outer radius of 18 km and 30 km respectively and 
a thickness of 5 km. The velocity was kept at 
v = 0.2c, but the inclination was decreased to 
25° . The fit shown in FigJTT] also reproduces the 
observed width of the line quite well. As with 
the spherical case, the "bump" at ~17 keV is well 
replicated in this model. 

We know that when a rotating disk is seen edge 
on, the line is split into a red and blue component. 
When relativistic beaming is significant, the red 
component is less intense than the blue. We next 
consider that the observed feature ~9 keV may 
in fact be this red component of the line split in 
the velocity field. If this were true, the actual line 
center would not be at 13.07 keV, but ~11 keV. 
The feature at ~9 keV would then be the red side 
of the split line, and the main feature at 13.07 keV 
would be the blue side of the feature. In order to 
reduce the number of free parameters, we assume 
the disk is rotating with Keplerian velocity, that 



v(r) 



GM 



(14) 



With these criteria, we now have four free pa- 
rameters for our model: inner and outer radius 
of the disk, mass of the compact object, and the 
inclination of the disk. Our best fit model with 
these parameters is a disk with inner and outer 
radius of 12 km and 16 km respectively, inclined 
at 32° around a compact object with M = 2M . 
The geometrical mesh representation of this model 
is shown in Fig |12| and the resulting spectrum in 
the bottom right panel of Fig{l3] convolved with 
a Gaussian of width 1.5 keV to match the resolu- 
tion of the observations. Not only is the main 
feature at 13.07 keV fit very well, but the two 
smaller features at ~ 9keV and ^17keV are as 
well. The feature at 9 keV is the red component 
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of the split line, the feature at 13.07 keV is the 
blue component, and the feature at ~17 keV is 
the contribution from the blue side of the weaker, 
but more abundant, Rb lines in this region. To 
show the breakdown in processes contributing to 
the line shape, we provide Fig |13[ a four panel fig- 
ure showing the spectrum with increasing levels of 
physical accuracy. We use 1000 frequency bands to 
show the fine structure of the spectrum. The top 
left panel shows the original non-shifted spectrum 
composed of several Rb transitions. The top right 
panel shows the effects of adding the gravitational 
field; the line is shifted to the red and slightly 
broadened. The bottom left is the addition of a 
Keplerian velocity profile, which has the effect of 
splitting the line into a red and blue component. 
Finally the bottom right panel adds the relativis- 
tic beaming effect which suppresses the red and 
enhances the blue component of the split line. We 
should note that the observed features at ~9 keV 
and ~17 keV are not necessarily real as they have 
low significance. Further, more sensitive observa- 
tions of XTE J1810-197 are therefore needed to 
obtain better data which may resolve this issue. 

These simulations have shown that atomic 
lines, specifically those of Rb, can explain the 13 
keV feature seen in XTE J1810-197. If the feature 
at ~9kEv and ^17keV are real, the model that 
best fits the data is a Keplerian disk with an inner 
and outer radius of 12 km and 16 km respectively 
in a gravitational field generated by a compact 
object of mass M = 2Mq. If these two smaller 
features are not real, the main feature at 13.07 
keV can be well reproduced by a rotating disk, 
rotating sphere or a static extended atmosphere 
all in a weaker gravitational field. 

5. Discussion 

We have investigated how the 13 keV feature 
of XTE J1810-197 may be reproduced under dif- 
ferent conditions and processes. We now take this 
knowledge and apply it to the possible theories 
for AXPs. For each emission process (cyclotron or 
atomic) we examine the feasibility for each AXP 
model. 

5.1. Cyclotron Emission 

Under the right conditions, both electron and 
proton cyclotron emission can reproduce the 13 



keV feature in XTE J1810-197. Within a dipole 
field, we found that emission from a thin region 
above the poles is sufficient. For proton cyclotron 
radiation, the B field within this region must be 
~2x 10 15 G and for electron cyclotron radiation, 
it must be - 1 x 10 12 G. 

In the Magnetar model, the surface is expected 
to have magnetic fields B ~ 1 X 10 14 -> 1 X 10 15 G. 
Cyclotron emission from the surface must there- 
fore be proton in nature in order to account for 
the line center of the 13 keV feature. This would 
support the theory that the burst originates from 
fractures in the neutron star crust. If the emission 
is in fact electron cyclotron emission, it must come 
from a region considerably higher than the neu- 



tron star surface. Fig 14 (blue Xs) shows the pos- 



sible electron cyclotron emission regions around 
a neutron star with a maximum field strength of 
1 x 10 15 G on its surface (~ 12 km). We can 
see that the possible emission regions are between 
~90 km and ~120 km. This would support the 
idea that the burst was caused by a reconnection 
event in the magnetosphere. It is difficult to de- 
termine whether the feature is caused by proton or 
electron cyclotron emission since the width of the 
line is dominated by the variability in the mag- 
netic field, and not by any differences between the 
two particles. The spin down rate of XTE J1810- 
197 implies a magnetic field of B ~ 3 x 10 14 G 
where the proton cyclotron interpretation requires 
a much larger field of -B ~ 2xl0 15 G. This discrep- 
ancy may therefore favour the electron cyclotron 
model. 

The magnitude of the magnetic field in the FBD 
model is expected to be < 1 x 10 13 G . This is 
much too low for proton cyclotron emission at the 
required energies and therefore only electron cy- 
clotron radiation from the surface of the neutron 
star is feasible in this model. Electron cyclotron 
radiation from the disk is ruled out since the field 
at the disk location (~10,000 km) is much too 
weak. 

In the QN model, only electron cyclotron emis- 
sion is plausible. The QN remnant is a quark star 



in a color-superconducting state ( Vogt et al. 2004 



Ouyed et al.||2005| a nd is not expected to have a 
crust ( jAlford et al.|[l999| . So although the field 
strengths at the surface (B ~ 1 x 10 15 G) are 
enough to provide the necessary proton cyclotron 
energies, no emission from this region will take 
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Electron cyclotron emission from the Kep- 

(red 



14 



place. 

lerian ring, however, may be possible. Fig 
Xs) shows the possible electron cyclotron emission 
regions around a quark star with a maximum mag- 
netic field of B ~ 1 X 10 15 G on its surface (radius 
of 8 km). If electron cyclotron emission around a 
quark star is responsible for the 13 keV feature, it 
must be emitted from a region between ~60 km 
and ~83 km. The Keplerian ring in the QN model 
is expected to be located between ~20 km and 
~100 km, which makes electron cyclotron radia- 
tion in this scenario a real possibility. 

Cyclotron emission in each model therefore ap- 
pears to be a viable cause of the 13 keV feature 
seen in XTE J1810-197. If XTE J1810-197 was the 
only available spectrum, this would be an appro- 
priate conclusion. However, emission lines have 
been observed in two other AXPs, 4U 0142+61 
and IE 1048.1, which also show a significant fea- 
ture at ~13 keV (Fig 15). As we have seen from 
our simulations, the location of the cyclotron line 
depends on the magnetic field strength within the 
emitting region and the inclination of that region 
with respect to the observer. For the feature to ap- 
pear in each AXP spectrum, the emission region in 
each must therefore have the same magnetic field 
and orientation. Or, which is even more unlikely, 
a suitable combination of the two that coinciden- 
tally results in the same line center. 

5.2. Atomic Line Emission 

We have shown that under various conditions 
(i.e. rotating disk, rotating shell, extended atmo- 
sphere) atomic line emission from Rubidium can 
adequately reproduce the 13 keV feature of XTE 
J1810-197. 

Rubidium 

Rubidium isotopes 85 Rb and 87 Rb are produced 
through the slow (s-) and rapid (r-) nuclear pro- 
cesses. The s-process occurs in low mass asymp- 



totic giant branch stars (AGB) (Garcia- Hernandez 



et al 2006) and the r-process likely during a su- 



pernova or quark nova explosion. In a supernova, 
much of the Rb is distributed into the interstellar 
medium, whereas in a quark nova the Rb will be 
retained in the Keplerian ring (depending on the 
quark nova parameters). 

In the case of a magnetar, most Rb is likely 



removed from the system, although some may re- 
main within the crust. If the "starquake" produces 
fractures in the neutron star crust, this Rb may be 
able to make its way to the surface and emit. The 
gravitational field produced by the neutron star 
must be such that the Rb line is shifted from its 
emitted energy of 14.299 keV to the observed en- 
ergy at 13.07 keV. The mass of the neutron star 
required for the Rb emission from a typical radius 
of 10 km is - O.6M . 

In the FBD model, the Rb expelled during the 
SN explosion may fall back and make up some of 
the disk. It is possible, therefore, for Rb to be 
present in the spectrum. For the accretion sce- 
nario, the shift of the line requires a neutron star 
with a mass of ~ 0.6M Q for an emitting radius 
(surface) of 10 km. 

The QN is an explosion of the neutron star crust 
resulting in an abundance of iron and neutrons. 
The r-process is therefore very efficient during a 
QN explosion. We have performed calculations 



using the code r-Java (Charignon et al. 2011) to 



calculate the abundance of Rb under the QN sce- 
nario and find that Rb is readily produced in the 
QN ejecta. Since the QN ejecta form the Keple- 
rian ring, we can expect a significant abundance 
of Rb available to emit. 

Broadening 

In order to explain the breadth of the 13 keV 
feature with Rb atomic lines, the atoms must be 
either moving a significant fraction of the speed of 
light, or extended around a compact massive ob- 
ject. For the Magnetar model, the Rb from pre- 
sumably the surface would have to be travelling 
~0.2c. Given that the rotation period of XTE 
J1810-197 is ^5.5 s, the radius of the neutron star 
would have to be ^5200 km to attain these speeds 
on its surface. 

This same argument holds for the accretion 
onto the neutron star in the FBD model. The 
emission, however, may come from the disk itself. 
If the disk were Keplerian, we can see from equa- 
tion [14] that for a velocity of 0.2c and an inner 
radius of ^10,000 km the mass of the compact 
object would have to be at least 196M Q , again 
unreasonable. 

The Rb emission in the QN model is expected 
from a Keplerian ring with an inner and outer ra- 
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dius between ~20 km and 100 km respectively. We 
have already modelled this situation, and found 
that a Keplerian disk with an inner and outer ra- 
dius of 12 km and 16 km respectively around a 
compact object with a mass M = 2Mq fits the 



observations very well (Fig 13). 

Gravitational broadening of the Rb lines re- 
quires an extended atmosphere. For a compact 
object with a mass of 1.4M Q , the inner and outer 
edge of the atmosphere must be ~ 13 and >^ 300 
km respectively. The emission, if the breadth is 
due to gravitational broadening, can therefore not 
come from the surface of the neutron star. This 
precludes the FBD theory, assuming the emission 
is from accretion onto the star. Emission from the 
disk (inner radius ~ 10,000 km) is also unlikely be- 
cause the mass of the compact object required to 
broaden the line would be ~ 714Af Q . In the Mag- 
netar model, emission of Rb from the surface of the 
star is discarded using the same argument. How- 
ever, emission may come from the atmosphere well 
above (~100 km) the magnetar surface. Ibrahim 
et al (2001) suggest that radiation pressure may 
lift a thin layer of material away from the mag- 
netar surface. Whether or not enough Rb could 
be driven off the surface, or if it could even reach 
an altitude of ^100 km or more is uncertain. The 
observation of the 13 keV feature in the tail of 
the burst does, however, seem to support the idea. 
Because of this we cannot fully rule out gravita- 
tional broadening in the Magnetar model. The QN 
model predicts an inner disk radius close to the re- 
quired value of ~13 km. Ouyed et al. ( 2007b ) pro- 



vides an expression for the outer radius of the Ke- 
plerian ring as a function of equilibrium tempera- 
ture and age. For an age of i « 8720 yrs and an 
equilibrium ring temperature of 0.5 keV, we calcu- 
late an outer radius of ~306 km which agrees with 
the required radius of >^300 km. Gravitational 
broadening in the QN model can therefore explain 
the full width of the 13 keV feature in XTE J1810- 
197, although it is more likely that both Doppler 
and Gravitational broadening contribute. 

6. Conclusion 
6.1. Summary 

We have performed an analysis of different ge- 
ometrical and physical conditions required to ex- 
plain both the location and shape of the 13 keV 



feature seen in XTE J1810-197. Both cyclotron 
(positron and electron) and Rb atomic line emis- 
sion are possible explanations for this particular 
feature. 

If the feature is proton cyclotron emission, the 
only model that can satisfactorily explain it is the 
Magnetar model. In this case the emission must be 
coming from the surface of the neutron star. Elec- 
tron cyclotron emission can be explained in the 
Magnetar, FBD, and QN models. In the Magnetar 
model, the electron emission must come from well 
above the neutron star, perhaps in the magneto- 
sphere. Due to the low surface magnetic fields ex- 
pected in the FBD model, electron cyclotron radi- 
ation would have to come from the stellar surface. 
In the QN model, the electron emission would 
come from the Keplerian ring where the magnetic 
fields are lower. 

Atomic line emission from Rb is also possible 
in each model. Rb is the most likely species re- 
sponsible for the 13 keV feature, and is expected 
to exist, at least minimally, in each theory. The 
Magnetar and FBD models predict small traces of 
Rb whereas the QN theory predicts large quanti- 
ties. The breadth of the Rb line requires velocities 
approaching the speed of light or extended atmo- 
spheres, both unexpected in the Magnetar model. 
The fall-back disk of the FBD model is too large 
for either of these processes to adequately explain 
the shape of the 13 keV feature. This does not, 
however, preclude the presence of a fall-back disk 
around this, or any other AXP source. The FBD 
model requires an additional mechanism to ex- 
plain the bursting, which may also provide a site 
for the emission of the 13 keV feature. The QN 
model, on the other hand, predicts a close, rapidly 
rotating Keplerian ring which is more than ade- 
quate to explain the breadth in terms of Doppler 
broadening alone. Gravitational broadening may 
also contribute to the breadth of the line in this 
situation. The only theory that can satisfactorily 
account for all properties of the 13 keV feature un- 
der the atomic line hypothesis is that of the quark 
star in the context of the QN explosion. 

If XTE J181 0-197 were the only source with 
the 13 keV feature, cyclotron emission would be 
a very likely explanation. However, in every AXP 
with an emission spectrum a feature with the same 
line center has been observed. This would require 
all three sources to have emission regions with the 
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same magnetic field and orientation with respect 
to earth. Although possible, this scenario is not 
very likely. For this reason we do not believe that 
the 13 keV feature in XTE J1810-197 is cyclotron 
in nature, and is more likely to be an atomic tran- 
sition of Rb. A summary of our finding is pre- 
sented in Tabic [2] 

6.2. Predictions 



We have found that the QN model provides the 
best explanation for the location and breadth of 
the 13 keV feature in XTE J1810-197. We can 
therefore make several predictions based on this 
model: 



field penetrates the ring atmosphere, we ex- 
pect to see electron cyclotron lines along- 
side atomic lines in the spectra of future 
AXP observations. 



• The width of the 13 keV feature i n XTE 
J181 0-197 is not well constrained ( |Woods| 



et al.| [2005). Future, high resolution, obser- 



vations will provide better data allowing us 
to more accurately determine the required 
velocities. This will put strong constraints 
on AXP models, and provide invaluable in- 
formation regarding the environment around 
these enigmatic objects. 



• The Keplerian ring in the QN model is 
unique to AXPs and will split the atomic 
lines into a red and blue component. SGRs, 
under this model, are explained by a co- 
rotating shell at much lower velocities. We 
therefore do not expect to see any line split- 
ting in SGRs. 

• We expect to see a feature at ~17 keV rep- 
resenting the less intense, but more abun- 
dant Rb lines. A strong detection of this 
feature would support this model. A null de- 
tection, however, may not rule it out since 
the strength of these lines depends on the 
excitation conditions (i.e. temperature and 
density). 

• The relative strength of the ~17 keV feature 
and the main ~13 keV feature is dependent 
on the relative abundance of RbXXXVI and 
RbXXXVII (the former has no lines around 
17 keV). A strong detection of the ~17 keV 
feature will therefore provide constraints on 
the electron density and temperature of the 
emitting region. 

• If atomic lines are responsible for the 13 keV 
feature, we should see evidence of other lines 
during quiescence, in the form of absorption 
features. 

• During and shortly after a burst, the atmo- 
sphere around the Keplerian ring is expected 
to be abundant in both ionized heavy- 
elements and electrons. Since the magnetic 
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Fig. 2. — Left: Cyclotron emission lines in a constant magnetic field of 1 x 10 12 G and 2.1 x 10 15 G for 
the electron and proton cases respectively. Angles at 30° and 60° between the magnetic field and observer 
are shown, as well as the observed spectrum of XTE J1810-197. Right: Proton cyclotron emission lines in 
a variable magnetic field; B — 2.1 x 10 15 G ±dB. Two values of dB are shown, as well as the observed 
spectrum of XTE J1810-197. It is clear that the larger variation in magnetic field produces the broader line. 




Fig. 3. — Cyclotron emission from a spherical shell within a dipolc magnetic field. Two different shell widths 
are plotted, 1km and 2km. The outer radius of the shell in both cases is 12 km. The observed spectrum of 
XTE J1810-197 is plotted in the background for reference. Left: Electron cyclotron model. Right: Proton 
cyclotron model. 
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Fig. 4. — Electron cyclotron emission from the polar regions between (f> = 0° and <j> — 20° of a spherical shell 
within a dipole magnetic field. Two different shell widths are plotted, 1 km and 2 km. The outer radius of 
the shell in both cases is 12 km. The observed spectrum of XTE J1810-197 is plotted in the background for 
reference. Left: Electron cyclotron model. Right: Proton cyclotron model. 




Fig. 5. — Electron (left panel) and proton (right panel) cyclotron emission from the polar regions between 
= 0° and = 20° of a spherical shell within a dipole magnetic field. Plotted are the results of simulations 
of four different angles between the magnetic moment of the dipole and the observer. In all cases the thin 
atmosphere of 1 km is used. The observed spectrum of XTE J1810-197 is plotted for reference. 
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Fig. 6. — Best fit model to the observed 13 keV feature of XTE J1810-197 for cyclotron emission. The 
emission is from a thin spherical shell with an inner and outer radius of 11.9 km and 12 km respectively, 
centred around the pole (4> = 0° — > 20°) of a magnetic dipole field. Left: Electron cyclotron model with a 
dipole magnetic moment of jjb m = 1 x 10 27 A m 2 at an inclination of 66° to the observer. Right: Proton 
cyclotron model with a dipole magnetic moment of // m = 2 x 10 30 A m 2 at an inclination of 56° to the 
observer. Each spectrum is convolved with a Gaussian of width 1.5 keV to simulate the energy resolution of 
the observations. 
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Fig. 8.— Model spectrum of Rb transitions at T = 1.5 x 10 s K. Rb XXXVII and Rb XXXVI transitions are 
shown. The abundance of RbXXXVI is 10% that of RbXXXVII. 1000 frequency bands were used to show 
the fine structure. Left: Original, unshifted spectrum. Right: Gravitationally shifted lines from a IAMq 
compact object. The emission is from a spherical shell with an inner and outer radius of 24.5 and 25.5 km 
respectively. 




Fig. 9. — Constant velocity, Rb atomic line model. Left: A spherical shell with an inner and outer radius of 
24.5 and 25.5 km respectively. Right: A disk with an inner and outer radius of 24.5 and 25.5 km respectively. 
The rotational velocity is v — 0.2c in both cases. Two different inclination of the rotational plane are plotted 
for each model. 
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Fig. 10. — Gravitational broadening from a disk in a gravitational field. The inner and outer radius of 
the disk is 10 km and 300 km respectively, and the mass of the compact object is 2M Q . a) Broadening of 
the line with a constant density profile throughout the disk, b) Broadening of the line with a 1/r 2 density 
profile, c) Broadening of the line with a 1/r 3 density profile, d) Broadening of the line with a 1/r 3 profile 
but convolved with a Gaussian of width 1.5 keV to simulate the energy resolution of the observations. 
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Fig. 11. — Best fit models under the atomic line model for a rotating spherical shell (red) and rotating disk. 
In each model, a constant rotational velocity of v — 0.2c is used. The data points were convolved with a 
Gaussian of width 1.5 keV to match observations. 




Fig. 12. — Geometrical representation of the Keplerian disk model. The mass of the central object is 2M Q , 
the inner and outer radius of the disk is 12 km and 16 km respectively. The system is inclined at 32°. Two 
different viewpoints are shown, the left panel from the front (observer) and the right panel from the side. 
The compact object is shown with a radius of 8 km, and the grid is spaced in 10 km intervals. 
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Fig. 13. — Sequence of panels showing the evolution of the line shape of the best fit atomic line as more 
physical processes are added. 1000 frequency bands were used to show the fine structure of the spectrum. 
The model consists of a Keplerian disk with an inner and outer radius of 12 and 16 km respectively, an 
inclination of 32° and a compact object of 2M . Top Left: Original, unmodified line shape contributed to 
by Rb transitions. Top Right: Gravitational field is added, causing the line to shift and broaden. Bottom 
Left: Keplerian velocity is added to the disk causing the line to split into a red and blue component of 
roughly equal intensity. Bottom right: Application of relativistic beaming causes the red side of the split 
line to be suppressed while the blue is enhanced. The final panel represents our best fit Keplerian disk model, 
convolved with a Gaussian of width 1.5 keV to simulate the energy resolution of the observations. 
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Fig. 14. — Plot of magnetic field strength with distance from an ideal magnetic dipole. Two different field 
strengths are shown. The blue xs represent the Magnetar model where the maximum magnitude at 12 
km (assumed neutron star surface) is 1 x 10 15 G. The red dots represents a quark star with a maximum 
magnitude of 1 x 10 15 G at its surface of 8 km. The green dotted line shows the required field for electron 
cyclotron emission to account for the position of the 13 keV feature. We can see that the emission region 
for the QN model is closer in than for the Magnetar model. 



25 




2G 



Table 1 

Atomic transitions between ~9 keV and ~18.5 keV from the NIST database. 



Species 


Energy fkeVl 


A 


Config 


Terms 






Rb XXXVII 


18.596 


2.24E+13 


Is 


- 7p 


2S - 


2P* 


2 - 


2 


Rb XXXVII 


18.457 


3.58E+13 


Is 


- 6p 


2S - 


2P* 


2 - 


2 


Rb XXXVII 


18.231 


6.40E+13 


Is 


- 5 P 


2S - 


2P* 


2 - 


4 


Rb XXXVII 


18.225 


6.27E+13 


Is 


- 5 P 


2S - 


2P* 


2 - 


2 


Rb XXXVII 


17.808 


1.27E+14 


Is 


- 4p 


2S - 


2P* 


2 - 


4 


Rb XXXVII 


17.797 


1.25E+14 


Is 


- 4p 


2S - 


2P* 


2 - 


2 


Rb XXXVII 


16.896 


3.11E+14 


Is 


- 3p 


2S - 


2P* 


2 - 


4 


Rb XXXVII 


16.870 


3.10E+14 


Is 


- 3 P 


2S - 


2P* 


2 - 


2 


Rb XXXVII 


14.299 


1.15E+15 


Is 


- 2p 


2S - 


2P* 


2 - 


4 


Rb XXXVII 


14.209 


1.18E+15 


Is 


- 2r> 


2S - 


2P* 


2 - 


2 


Rb XXXVI 


13.891 


1.66E+15 


ls2 - 


ls.2p 


IS - 


IP* 


1 - 


3 


Rb XXXVI 


13.794 


4.58E+14 


ls2 - 


ls.2p 


IS - 


3P* 


1 - 


3 


Ni XXVIII 


10.356 


2.18E+13 


Is 


- 5 P 


2S - 


2P* 


2 - 


4 


Ni XXVIII 


10.115 


4.33E+13 


Is 


-4p 


2S - 


2P* 


2 - 


4 


Ni XXVII 


9.899 


3.35E+13 


ls2 - 


ls.5p 


IS - 


IP* 


1 - 


3 


Ni XXVII 


9.897 


5.20E+12 


ls2 - 


ls.5p 


IS - 


3P* 


1 - 


3 


Ni XXVII 


9.675 


6.38E+13 


ls2 - 


ls.4p 


IS - 


IP* 


1 - 


3 


Ni XXVII 


9.670 


1.00E+13 


ls2 - 


ls.4p 


IS - 


3P* 


1 - 


3 


Co XXVII 


9.622 


1.88E+13 


Is 


- 5p 


2S - 


2P* 


2 - 


4 


Ni XXVIII 


9.594 


1.06E+14 


Is 


-3p 


2S - 


2P* 


2 - 


4 


Ni XXVIII 


9.585 


1.06E+14 


Is 


-3p 


2S - 


2P* 


2 - 


2 


Co XXVII 


9.398 


3.73E+13 


Is 


-4p 


2S - 


2P* 


2 - 


4 


Ni XXVII 


9.191 


1.63E+14 


ls2 - 


ls.3p 


IS - 


IP* 


1 - 


3 


Co XXVI 


9.182 


2.92E+13 


ls2 - 


ls.5p 


IS - 


IP* 


1 - 


3 


Co XXVI 


9.180 


4.00E+12 


ls2 - 


ls.5p 


IS - 


3P* 


1 - 


3 


Ni XXVII 


9.180 


2.40E+13 


ls2 - 


ls.3p 


IS - 


3P* 


1 - 


3 


Co XXVI 


8.975 


5.68E+13 


ls2 - 


ls.4p 


IS - 


IP* 


1 - 


3 


Co XXVI 


8.971 


7.70E+12 


ls2 - 


ls.4p 


IS - 


3P* 


1 - 


3 


Fe XXVI 


8.916 


1.61E+13 


Is 


- 5p 


2S - 


2P* 


2 - 


4 


Co XXVII 


8.914 


9.17E+13 


Is 


-3p 


2S - 


2P* 


2 - 


4 
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Table 2 

Summary of possible explanations for the 13 keV feature of XTE J1810-197 and how each 

AXP model performs. 



Model 


Cyclotron Emission* 




Atomic Line Emission 




Proton 


Electron 


Rubidium 


Doppler Broadening 


Gravitational Broadening 


Magnet ar 


Yes 


Yes 


Unlikely 


No 


Unlikely 


Fall-Back Disk 


No 


Yes 


Unlikely 


No 


No 


Quark Nova 


No 


Yes 


Yes 


Yes 


Yes 



*The feasibility of each model under cyclotron emission is for the single case of XTE J1810-197, and 
does not take into account the argument of multiple sources with the same feature. 
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